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ABSTRACT 

From a set of stellar spectropolarimetric observations, we report the detection of sur- 
face magnetic fields in a sample of four solar-type stars, namely HD 73350, HD 76151, 
HD 146233 (18 Sco) and HD 190771. Assuming that the observed variability of polari- 
metric signal is controlled by stellar rotation, we establish the rotation periods of our 
targets, with values ranging from 8.8 d (for HD 190771) to 22.7 d (for HD 146233). 
Apart from rotation, fundamental parameters of the selected objects are very close to 
the Sun's, making this sample a practical basis to investigate the specific impact of 
rotation on magnetic properties of Sun-like stars. 

We reconstruct the large-scale magnetic geometry of the targets as a low-order 
(JL < 10) spherical harmonics expansion of the surface magnetic field. From the set 
of magnetic maps, we draw two main conclusions, (a) The magnetic energy of the 
large-scale field increases with rotation rate. The increase of chromospheric emission 
with the mean magnetic field is flatter than observed in the Sun. Since the chromo- 
spheric flux is also sensitive to magnetic elements smaller than those contributing to 
the polarimetric signal, this observation suggests that a larger fraction of the sur- 
face magnetic energy is stored in large scales as rotation increases, (b) Whereas the 
magnetic field is mostly poloidal for low rotation rates, more rapid rotators host a 
large-scale toroidal component in their surface field. From our observations, we infer 
that a rotation period lower than ^12 days is necessary for the toroidal magnetic 
energy to dominate over the poloidal component. 

Key words: Line : polarization - Stars : rotation - activity - magnetic fields - Star : 
individual : HD 73350 - HD 76151 - HD 146233 - HD 190771 



1 INTRODUCTION 

According to dynamo models, the variable magnetic field 
of the Sun is the consequence of the interplay between two 
main ingredients. The first ingredient is the vertical and lat- 
itudinal differential rotation that succeeds at generating a 
large-scale toroidal magnetic field from an initial poloidal 
field. The second ingredient is still a matter of debate, with 
models invoking either the cyclonic convection in the con- 
vection zone (Parker 1955) or the transport of decaying ac- 
tive regions by meridional circulation (Dikpati et al. 2004) 
as possible processes to regenerate the poloidal magnetic 
component. When acting together, both effects succeed at 
building continuously a large-scale magnetic field that oscil- 



lates with time, giving rise to the 22 yr period of the solar 
cycle. Despite considerable progress in this field since the 
very first solar dynamo model proposed by Parker (1955), 
there are still many aspects of solar magnetism that the cur- 
rent models cannot reproduce or did not thoroughly explore 
(see, e.g., the reviews of Ossendrijver 2003 or Brandenburg 
& Subramanian 2005). 

Our understanding of the solar dynamo can benefit from 
the observation of solar-type stars, where dynamo types 
marginal or inactive in the Sun can be observed, either be- 
cause these analogues of the Sun are caught by chance in an 
unfrequented activity state or because their physical proper- 
ties (in particular their mass and rotation rate) differ suffi- 
ciently from the Sun's to lead to a different dynamo output. 
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Chromospheric emission, observed on stars possessing an ex- 
tended convective envelope, has been monitored for decades 
as an indirect tracer of stellar magnetism (e.g. Baliunas et 
al., 1995) and offered a first opportunity to investigate the 
solar-stellar connexion. A variety of behaviours are observed, 
from erratic variations (common in young solar analogues) 
to smooth cycles (which seem typical of older stars like the 
Sun). Yet, the total chromospheric flux provides little in- 
formation about the spatial organization of the magnetic 
field at the stellar surface, the key observable to collect for a 
comparison with magnetic geometries predicted by dynamo 
models. 

Any useful effort at drawing a solar- stellar connexion 
from direct magnetic field measurements needs to address 
the case of solar twins, i.e. objects with stellar parameters 
very close to the Sun's, including a similar level of magnetic 
activity. Despite the strong scientific motivation, measuring 
a magnetic field on a strict twin of our Sun has always failed 
so far. Although being successful on a few objects quite dif- 
ferent from the Sun by their mass (Valenti et al. 1996, Rein- 
ers &; Basri 2007) or rotation rate (Riiedi et al. 1997, Donati 
et al. 2003, Petit et al. 2005), measurements of Zeeman sig- 
natures have always failed at detecting the magnetic field of 
low-activity solar twins, for well-identified reasons. 

A first option at our disposal, consisting in measuring 
the magnetic broadening of spectral lines (Saar 1988) is sen- 
sitive to the total photospheric magnetic flux but is plagued 
by the small fraction of the photosphere covered by strong 
magnetic fields. A second option consists in analyzing circu- 
larly polarized signatures inside line profiles (Donati et al. 
1997) to recover information on the line-of-sight component 
of the magnetic vector. This second strategy suffers from the 
complex magnetic topology of cool stars, featuring a mixture 
of opposite polarities on the visible hemisphere, so that their 
respective polarized signatures cancel out whenever the ro- 
tational broadening is not sufficiently large to disentangle 
their respective contributions. For slow-rotators, only the 
largest-scale components of the surface field can add-up con- 
structively to produce detectable circular polarization. On 
the Sun, this global component of the magnetic field dis- 
plays a strength limited to a few Gauss only (Babcock & 
Babcock 1955, Smith & Balogh 1995, Sanderson et al. 2003, 
Mancuso & Garzelli 2007). On other low-activity stars, the 
very tiny Zeeman signatures associated to such global fields 
have always escaped the scrutiny of observers so far. 

We achieve this detection in the present study, with 
the help of spectropolarimetric data sets collected with the 
NARVAL spectropolarimeter. From a set of observations of 
a sample of four Sun-like stars covering a range of rotation 
periods, we reconstruct their large-scale photospheric mag- 
netic geometry and discuss the impact of rotation on their 
magnetic properties. 



2 STELLAR SAMPLE 

Our stellar sample is constituted of four nearby dwarfs. 
Their physical parameters are chosen to be as close as pos- 
sible to the Sun's, except for their rotation period (Table 
[TJ . According to the spectral classification study of Valenti 
& Fischer (2005), all selected objects have an effective tem- 
perature compatible with the solar value within 1.3cr. They 



also possess a surface gravity compatible with the Sun, ex- 
cept HD 76151 (lying 1.8cr above the solar value). Two of 
our targets are significantly over-metallic compared to the 
Sun, with HD 190771 differing from the Sun by as much as 
4.5a. 

Evolutionary models matching these atmospheric pa- 
rameters provide us with estimates of stellar ages and masses 
(also listed in Valenti & Fischer, 2005). Stellar masses are 
all close to solar, except HD 76151 that lies 2a above 1 M©. 
Note however that Nordstrom et al. (2004) propose for this 
star a mass of 0.92l°;fg M©, while all other stars of our list 
have consistent mass estimates in both catalogues. All ages 
are consistent with a solar age, but the related uncertainties 
are generally in excess of 2 Gyr. In particular, HD 190771 
is probably younger than the Sun (given its metallicity and 
fast rotation), but is still sufficiently evolved to be taken as 
a probable main-sequence dwarf. As far as the dynamo itself 
is concerned, this large uncertainty in stellar ages could be 
a relevant issue in case the depth of the convection zone was 
expected to significantly vary during the course of the main 
sequence. 1M© models produced by stellar evolution codes 
predict however that the radius of the base of the convection 
zone remains constant within 5% over the first 9 Gyrs of the 
main sequence (S. Theado, private communication). 

From this set of stellar parameters, and given the re- 
lated uncertainties, we conclude that the internal structure 
of the selected objects are probably very close to the Sun's, 
with a dynamo action that should differ from the solar dy- 
namo through the dominant effect of rotation, including dif- 
ferential rotation (and possibly an unknown meridional cir- 
culation) . 



3 SPECTROPOLARIMETRIC DATA SET AND 
MAGNETIC MAPPING 

The observational material consists of high-resolution spec- 
tra obtained simultaneously in classical spectroscopy (Stokes 
I) and circularly polarized light (Stokes V) in 2007 win- 
ter and summer, using the newly installed NARVAL stellar 
spectropolarimeter at Telescope Bernard Lyot (Observatoire 
du Pic du Midi, France). As a strict copy of ESPaDOnS (Pe- 
tit et al. 2003), NARVAL provides full coverage of the op- 
tical domain (370 nm to 1,000 nm) in a single exposure, at 
a resolving power of 65,000, with a peak efficiency of about 
15% (telescope and detector included). It is constituted of a 
bench-mounted spectrograph (based on a dual-pupil design 
and stored in a double thermal enclosure), fiber- fed from a 
Cassegrain-mounted module containing all polarimetric fa- 
cilities. A series of 3 Fresnel rhombs (two half- wave rhombs 
that can rotate about the optical axis and one quarter-wave 
rhomb) are employed to perform a very efficient polarimetric 
analysis over the whole spectral domain. They are followed 
by a Wollaston prism which splits the incident light into two 
beams, respectively containing light linearly polarized per- 
pendicular/parallel to the axis of the prism. The two beams 
produced by the Wollaston prism are imaged onto the two 
optical fibers that carry the light to the spectrograph. Each 
Stokes V spectrum is obtained from a combination of four 
sub-exposures taken with the half-wave rhombs oriented at 
different azimuths (Semel et al. 1993). The data reduction 
is performed by Libre- Esprit, a dedicated, fully automated 
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Table 1. Fundamental parameters of the stellar sample. The effective temperature (T e ff ), surface gravity log(g), metallicity [M/H], mass, 
age and projected rotational velocity (v sinz) are taken from Valenti & Fischer (2005), except the -usini value of HD 146233 (see Sect. 
3). The equatorial rotation period (P^ot)' difference in rotation rate between pole and equator (dfi) and inclination angle are derived 
from Zeeman-Doppler imaging. 



Name 


T e ff 

(K) 


log(g) 
[cm.s -2 ] 


[M/H] 
[Sun] 


Mass 
(M ) 


Age 
(Gyr) 


v s'mi 
(kms- 1 ) 


■peq 
^ rot 

(d) 


dQ 
(rad.d- 1 ) 


inclination 

(°) 


Sun 
HD 146233 
HD 76151 
HD 73350 
HD 190771 


5770 
5791 ± 50 
5790 ± 50 
5802 ± 50 
5834 ± 50 


4.44 
4.41 ±0.06 
4.55 ±0.06 
4.48 ± 0.06 
4.44 ±0.06 


0.00 
0.03 ±0.03 
0.07 ±0.03 
0.04 ±0.03 
0.14 ±0.03 


1.0 

0.98 ±0.13 
1.24 ±0.12 
1.01 ±0.14 
0.96 ±0.13 


4.3 ± 1.7 

4.7i 2 2 ; 7 

3.6iJ;S 
4 i+2-O 

^• 1 -2.7 
o 7+1-9 
z - ' -2.0 


1.7 

2.1 ±0.5 

1.2 ±0.5 
4.0 ±0.5 

4.3 ±0.5 


25 

22.7 ±0.5 
20.5 ±0.3 
12.3 ±0.1 
8.8 ±0.1 


0.05 

0.2 ±0.2 
0.12 ±0.03 


30 ± 15 

'°-20 

50 ± 10 
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Figure 1. From left to right, time-series of Stokes V profiles of HD 146233, HD 76151, HD 73350 and HD 190771, in the rest-frame 
of the star. Red lines represent the data and black lines correspond to synthetic profiles of our magnetic model. Successive profiles are 
shifted vertically for display clarity and rotational phases of observations are indicated in the right part of the plot. 



software described by Donati et al. (1997) and implement- 
ing the optimal spectral extraction principle of Home (1986) 
and Marsh (1989). 

A total of 9 to 14 spectra was collected for each star, at a 
rate of one spectrum every clear night (Table[2]). A single, av- 
erage photospheric line profile was extracted from each spec- 
trum using the LSD technique (Donati et al. 1997), accord- 
ing to a line-list matching a solar photospheric model. Using 
a total of some 5,000 atomic spectral lines with wavelengths 
between 370 nm and 1 /xm, the noise level of the mean Stokes 
V profiles is reduced by a factor of about 40 with respect to 
the initial spectrum. The resulting noise level, listed in Table 
[2] ranges between 2 x 10~ 5 and 8 x 10 ~ 5 I c (where I c de- 
notes the continuum level) depending on the source, on the 
fluctuating weather and on the adopted exposure time. The 
recorded sets of Stokes V profiles are plotted in Fig. [I] Cir- 



cularly polarized signals are detected for each star, centered 
at the radial velocity of the unpolarized mean line profile. 
We interpret this spectral line polarization as Zeeman sig- 
natures related to the presence of large-scale photospheric 
magnetic fields. This detection, achieved for the first time 
in solar twins, is made possible by the combination of the 
large spectral domain (ensuring highly efficient multi-line 
techniques) and high S/N of polarized spectra (achievable 
thanks to the high collecting power of the instrumental de- 
vice) . 

Assuming that the observed temporal variability of 
Stokes V profiles is controlled by the stellar rotation, we re- 
construct the magnetic geometry of our targets by means of 
Zeeman-Doppler imaging (ZDI). We employ here the mod- 
elling strategy of Donati & Brown (1997), including also the 
spherical harmonics expansion of the surface magnetic field 
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implemented by Donati et al. (2006) in order to easily dis- 
tinguish between the poloidal and toroidal components of 
the reconstructed magnetic field distribution. We use a lin- 
ear limb- darkening coefficient equal to 0.75 (Neckel 2003, 
but note that the imaging procedure is mostly unsensitive 
to small variations of this parameter) and we adopt the 
i>sini values derived by Valenti & Fischer (2005), except 
for HD 146233 for which we choose i>sinz=2.1 kms -1 (since 



Table 2. Journal of observations. From left to right, we list the 
Julian date, the exposure time, the error bar in Stokes V LSD 
profiles and the phase of the rotational cycle at which the obser- 
vation was made (calculated according to the rotation periods of 
Table 1 and imposing JD=2454101.5 as phase zero for all stars). 



Name 


Julian Date 


exp. time 


0LSD 


rot 


phase 




(2,450,000+) 


sec. 


io- 5 / c 






HD 146233 


4307.41 


3600.0 


5.894 





0711 




4309.36 


3600.0 


2.201 





1571 




4311.38 


3600.0 


2.029 





2459 




4312.37 


3600.0 


1.951 





2894 




4313.37 


3600.0 


3.395 





3336 




4315.37 


3600.0 


2.294 





4216 




4316.37 


3600.0 


2.174 





4658 




4322.37 


3600.0 


2.850 





7298 




4323.37 


3600.0 


3.234 





7740 




4331.35 


3600.0 


2.465 


1 


1258 


HD 76151 


4122.51 


2000.0 


7.982 





0250 




4126.55 


3200.0 


4.887 





2223 




4127.56 


3200.0 


4.178 





2712 




4128.57 


3200.0 


3.727 





3206 




4130.56 


3200.0 


2.956 





4176 




4132.55 


3200.0 


3.457 





5149 




4133.59 


3200.0 


3.570 





5655 




4134.56 


3200.0 


3.335 





6129 




4135.58 


3200.0 


4.114 





6626 




4136.55 


3200.0 


4.404 





7100 




4146.58 


3200.0 


3.954 


1 


1993 




4150.56 


3200.0 


3.471 


1 


3933 




4159.46 


3200.0 


3.573 


1 


8277 


HD 73350 


4127.47 


2400.0 


7.113 





1122 




4128.43 


2400.0 


6.274 





1902 




4130.43 


2400.0 


5.329 





3518 




4133.45 


2400.0 


8.391 





5983 




4134.47 


2400.0 


6.899 





6813 




4135.45 


2400.0 


4.943 





7605 




4136.42 


2400.0 


5.519 





8392 




4138.42 


3600.0 


5.659 


1 


0021 




4140.43 


2400.0 


5.569 


1 


1656 


HD 190771 


4307.52 


2400.0 


4.584 





4121 




4308.47 


2400.0 


4.313 





5194 




4309.50 


2400.0 


3.506 





6369 




4310.52 


2400.0 


3.159 





7528 




4311.52 


1600.0 


4.230 





8663 




4312.50 


1600.0 


4.055 





9783 




4313.50 


1600.0 


5.451 


1 


0909 




4315.51 


1600.0 


4.439 


1 


3199 




4316.51 


1600.0 


4.221 


1 


4332 




4321.48 


1600.0 


4.772 


1 


9983 




4322.50 


1600.0 


5.270 


2 


1146 




4323.50 


1600.0 


5.792 


2 


2281 




4327.50 


1600.0 


4.615 


2 


6830 




4331.43 


1600.0 


4.116 


3 


1291 




15 20 

Rotation period (d) 



Figure 2. Periodograms obtained for HD 73350, HD 76151 and 
HD 146233 (red, green and blue lines, respectively). 
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fi eq (rad/d) 

Figure 3. x 2 map obtained when varying the differential rotation 
parameters of HD 190771. f2 eq is the rotation rate of the equator 
and dQ represents the difference in rotation rate between pole 
and equator (assuming a sin 2 dependence of rotation rate with 
colatitude). 



the estimate of Valenti & Fischer is not compatible with the 
rotation period we derive for this star). Prior to rotational 
broadening, we adopt a very simple gaussian model for the 
unpolarized line profile, in which any type of line broad- 
ening other than rotational enter, including instrumental 
smearing. We only have to slightly adjust the FWHM of the 
local profile (between 10.1 and 10.4 kms -1 , depending on 
the star) to obtain a convincing adjustment of the observed 
Stokes I LSD profiles. 

To determine the rotation period of our targets, we use 
the principle of maximum entropy image reconstruction and 
calculate a set of magnetic maps for each star, assuming var- 
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ious values of the rotation period (Fig.JSJ. We then adopt the 
rotation period that minimizes the x of the reconstructed 
spectra, at fixed information content (following the approach 
of Petit et al. 2002). For all stars except HD 190771, the tem- 
poral evolution of the polarized signal is consistent with a 
solid-body rotation. For HD 190771 however, this assump- 
tion does not yield a convincing adjustment of the Stokes 
V time-series (% 2 = 1.5). A better fit to the profiles (% 2 = 
1.1) is obtained by including a solar- like differential rotation 
law in the imaging procedure (Petit et al., 2002 and Fig.|3|. 
Using the same method with HD 73350, we also obtain a x 2 
minimum in the Q e q-dQ plane, but with error bars too large 
to exclude a solid-body rotation. The two other stars do not 
show any % 2 minimum along the dQ, axis. 

For HD 73350, HD 76151 and HD 190771, our esti- 
mates of the rotation period (see Table [I]) are 20% to 25% 
larger than literature values derived from the less accurate 
period determination using chromospheric activity indica- 
tors (Noyes et al. 1984, Wright et al. 2004). For HD 146233, 
the level of chromospheric emission and the periodic mod- 
ulation of Zeeman signatures provide consistent estimates. 
The observed rotation phases are listed in Table [2] using 
the rotation periods of Table [l] and imposing JD=2454101.5 
(2007 Jan. 01 at 00:00 UT) as phase zero for all stars. 

From the combined knowledge of the rotation period 
and i>sini, we can finally deduce the inclination angle i of 
the rotation axis with respect to the line of sight (Table [TJ, 
assuming that for all targets the stellar radius can be taken 
equal to the solar radius (luminosities and effective tempera- 
tures taken from Valenti & Fischer, 2005, ensure that stellar 
radii of our targets are equal to the Sun's within «5%). Us- 
ing this set of stellar parameters, we then reconstruct pho- 
tospheric magnetic maps of our stellar sample (Figs. [4] and 
[5]). By doing so, all data sets are fitted down to the noise 
level (0.9 < x 2 < 1.1, depending on the star). We note that 
the hypothesis of a purely potential field is actually suffi- 
cient to adjust the Stokes V profiles of HD 146233 and HD 
76151 at the same % 2 5 without affecting the magnetic map. 
For HD 73350 however, assuming a purely potential mag- 
netic field does not allow the inversion procedure to reach 
better than x 2 = 1.5. The situation is even worse with HD 
190771, for which the final x 2 cannot be better than 5 with- 
out the addition of a toroidal magnetic field. We check that 
the magnetic maps remain essentially unaffected by the £ 
limit, provided that the magnetic expansion respects £ > 4 
for HD 73350, £ > 3 for HD 190771, £ > 2 for HD 146233 and 
HD 76151. For all four stars, we choose to limit the magnetic 
model to a spherical harmonics expansion of degree £ < 10. 

From the estimated spherical harmonics coefficients, we 
can derive a set of numerical quantities to characterize the 
reconstructed magnetic topologies (Table[3}. To evaluate the 
uncertainties on these numerical estimates, we vary the val- 
ues of stellar parameters (rotation period, vsini, inclination 
angle) over the width of the error bar on each individual 
parameter and reconstruct a magnetic map for each new 
combination of the parameters. The observed variations in 
the output magnetic quantities are quoted as error bars in 
Table [3] and give an estimate of the uncertainties on the 
reconstructed magnetic geometries. 



4 DISCUSSION 

4.1 Magnetic energy of the large-scale field 

The mean (unsigned) magnetic field of the reconstructed 
maps increases with the rotation rate (Figs. [4] [5] and [6] and 
Table [3}. This trend is in qualitative agreement with previ- 
ous studies investigating the rotational dependence of var- 
ious magnetic tracers, in particular the chromospheric flux 
(e.g. Noyes et al. 1984, Wright et al. 2004, Hall et al. 2007a) 
or the Zeeman broadening of magnetically-sensitive spectral 
lines (Saar 2001). From our set of spectropolarimetric obser- 
vations, we derive that the magnetic energy of the large-scale 
field increases by a factor of ~ 120 from HD 146233 (P^ot 
= 22.7 d) to HD 190771 (P^ t = 8.8 d). This evolution af- 
fects both poloidal and toroidal magnetic components, but 
the poloidal magnetic energy grows by a factor of about 40, 
while the toroidal magnetic energy is increased by a factor 
oflO 4 . 

From our sets of spectra, we can estimate the chro- 
mospheric emission of the selected stars from the R^k pa- 
rameter (defined as the ratio of the chromospheric emission 
in Call H & K line cores to the stellar bolometric emis- 
sion). Plotted against the mean magnetic field (Fig. |6]), the 
chromospheric emission is following a power-law of the form 
\og(R' UK ) = 0.32 log(B m ean) — 4.98. We obtain a very sim- 
ilar result when using literature values of Ruk instead of 
NARVAL values (R^k ^ #mean)- This dependence differs 
from solar observations, from which an increase like 5j; e 6 an 
is derived in the quiet Sun and active regions (Schrijver et 
al. 1989, Ortiz & Rast 2005). 

In the case of low vsini values, Zeeman-Doppler mag- 
netic maps are only sensitive to low-order field components, 
while the measured chromospheric flux includes the con- 
tribution of smaller-scale elements as well (active regions 
and plages), irrespective of the rotation rate. Since the spa- 
tial resolution achieved through ZDI increases with vsinz, 
one could then argue that higher-order magnetic compo- 
nents might be resolved in the most rapidly-rotating stars 
of our sample and might thus be included in the magnetic 
energy we measure, explaining why the mean magnetic field 
seems to grow much faster with rotation than the chromo- 
spheric flux. We actually observe that spherical harmonics 
modes with £ < 3 (spatially resolved for all stars of our 
sample) always contain more than 94% of the total pho- 
tospheric magnetic energy, even for stars with the largest 
i>sinz, and despite the fact that our magnetic modelling 
allows for higher-order magnetic components to be recon- 
structed up to £ — 10. In fact, the v sini range of our sample 
(from 1.2 to 4.3 kms -1 ) is sufficiently limited to ensure a 
consistent behaviour of the imaging procedure for all targets 
(given the spectral resolution of NARVAL, limited to 65,000 
in polarimetric mode), and therefore a consistent estimate 
of the magnetic energy of the large-scale field. 

From these arguments, we propose that the observed 
discrepancy with the solar case is genuine. A first possibility 
is that horizontal magnetic fields (in particular the promi- 
nant toroidal field components seen in Fig. [5]) are not able to 
penetrate the chromosphere and contribute to produce en- 
hanced Call emission. To test this idea, we have calculated 
the power law we obtain by considering the radial field com- 
ponent only (Fig. [6] blue symbols). By doing so, we only 
slightly change the exponent (with R^ oc B^ean) an d re- 
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Figure 4. Magnetic maps of HD 146233 and HD 76151 (left and right panel, respectively). Each chart illustrates the field projection 
onto one axis of the spherical coordinate frame. The magnetic field strength is expressed in Gauss. Vertical ticks above the charts indicate 
the observed rotational phases. Note that color scales are not the same for every star. 



main still very far from the steeper solar dependence. A 
second possibility is that a higher coverage of cool spots in 
fast-rotating Suns may contribute to alter the exponent in 
the power- law (which is valid for the quiet Sun and plages, 
not for sunspots). However, the dominant contributor to the 
Zeeman signatures we observe is probably not originating 
from cool spots (except possibly the spot penumbra), since 
the contrast with the rest of the photosphere will mostly 
hide the signal formed in these dark regions. By constrain- 
ing the temperature of magnetic regions observed in cool 
active stars, previous studies (Riiedi et al. 1997, Petit et 
al. 2005) confirm that the observed magnetic elements are 
hotter than the surrounding (non magnetic) material. Pe- 
tit et al. (2005) also show that Stokes V Zeeman signatures 
formed by the large-scale field of solar-type stars are asym- 
metric, with a difference of area and amplitude between the 
two lobes of the profile typical of solar faculae. 

We therefore propose a third interpretation and sug- 
gest that a larger fraction of the surface magnetic energy 
is stored in large scales as rotation grows, so that a larger 
fraction of the actual magnetic energy becomes measurable 
through ZDI, even when the field reconstruction is limited to 
large-scale structures. This conclusion is also supported by 
the observed increase of the mean magnetic field, for which 
a power- law of the form B mean = 2.8 x 10 3 P~ t 1-8 is derived 
from our observations. This increase is steeper than that 
derived from Zeeman broadening of spectral lines (a mea- 
surement technique also sensitive to smaller-scale magnetic 



elements, see Sect. [T]), from which Saar (2001) obtains an 
exponent equal to -1.2. 

Finally, we note that the observed increase of the mean 
magnetic field with the rotation rate is also in agreement 
with recent 3-D MHD simulations of fast-rotating Suns 
(Brown et al. 2007a, 2007b). In these simulations, the global 
field becomes stronger within the bulk of the convection zone 
when the rotation rate is increasing. This behaviour con- 
trasts with classical solar MHD models where weak mag- 
netic fields with little global-scale structure are produced 
and where fluctuating parts dominate near the surface (Brun 
et al. 2004). In other words, a way to organize and amplify 
magnetic fields in these global MHD simulations would con- 
sist in increasing the star rotation, what we also observe here 
with these four Sun-like stars. 

4.2 Large-scale toroidal component 

The second noticeable effect of rotation is to increase the 
fraction of the magnetic energy stored into a large-scale 
toroidal component of the surface magnetic field (Figs.|4]and 
[5]). While the toroidal component is negligible for HD 146233 
and HD 76151 (representing respectively 1% and 7% of their 
large-scale magnetic energy, Table [3]), the magnetic field be- 
comes predominantly toroidal for stars rotating faster (up 
to 65% of the magnetic energy of HD 190771). The absence 
of any significant toroidal component for slow rotators is 
unlikely to come from an intrinsic limitation of ZDI in the 
low vsini regime (where the sharp line profiles may favour 
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Figure 5. Same as Fig.fl] for HD 73350 (left panel) and HD 190771 (right panel). 



Zeeman signatures visible in the core of the line, i.e. ra- 
dial/meridional magnetic fields against the azimuthal field 
component). The main reason, already invoked in sect. 4.1 



is that the vs'mi range of our sample is sufficiently limited 
to ensure a consistent response of the imaging procedure 
over this domain. Another confirmation of this idea comes 
from observations of £ Bootis A by Petit et al. (2005), where 
a predominant toroidal field is reported despite a projected 
rotational velocity as low as 3 kms -1 (close to the vsini 
estimate proposed by Valenti & Fischer for HD 146233). 

The mainly poloidal geometries of HD 146233 and 
HD 76151 are reminiscent of the large-scale solar magnetic 
geometry and might be typical of low- activity dwarfs, at 
least for stellar masses close to the Sun's. For higher spin 
rates, the development of a toroidal field is reminiscent of 
previous ZDI studies of stars displaying a much higher mag- 
netic activity than the Sun (e.g. Donati et al. 2003, Petit 
et al 2005). Again, these observational results are now sup- 
ported by recent numerical simulations of dynamo action in 
rapidly rotating solar- type stars (Brown et al. 2007b), where 
an increase by a factor 3 of the rotation rate compared to the 
solar case results in a stable toroidal magnetic field at global 
scales, throughout the bulk of the convective envelope and 
without the help of any tachocline. The ordering of the mag- 
netic field onto a large-scale toroidal component when the 
star rotation increases is usually associated in MHD sim- 
ulations with the inclusion of a tachocline at the base of 
the convective zone (i.e. the addition of a strong local shear 
between the differentially rotating convective zone and the 
solid-rotating radiative core) . Solar models with a tachocline 



lead to the building of strong large-scale toroidal fields con- 
fined in this interface layer (Browning et al., 2006 & 2007), 
but the ability of these toroidal fields to travel the whole con- 
vective zone while keeping their identity remains neverthe- 
less a matter of debate (e.g., Dorch & Nordlund 2001, Jouve 
& Brun 2007) and the last simulations of rapidly rotating 
Suns without tachocline by Brown et al. precisely overcome 
this problem. From our analysis, and despite the small num- 
ber of objects in this first available stellar sample, we infer 
that a rotation period smaller than « 12 d is necessary for 
the toroidal component to reach an energetic level similar to 
that of the large-scale poloidal field. This is again in agree- 
ment with Brown et al. (2007b), where a dynamo model at 
Q = 3^0 produces a toroidal/poloidal field strength ratio 
greater than unity. 



4.3 Energy splitting between low-order magnetic 
components 

No clear rotational dependence is observed for the geomet- 
rical complexity of the reconstructed surface magnetic field. 
If we first consider the poloidal field only, a majority of 
the magnetic energy of this component is showing up in 
the dipole for HD 76151 and HD 190771 (with 79% and 
43% of the total poloidal magnetic energy, respectively) . For 
HD 146233, we reconstruct a predominantly quadrupolar 
field (56% of the poloidal magnetic energy). For HD 73350, 
the poloidal magnetic energy is more evenly distributed 
within the dipole, quadrupole and octopole (with 24%, 29% 
and 33% respectively). 
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If we now consider the toroidal field alone (whenever 
its contribution is not negligible, i.e. for HD 73350 and HD 
190771), we always measure a predominance of the mode 
(Z = 1, ra = 0) (87% and 81% of the total toroidal mag- 
netic energy of HD 73350 and HD 190771, respectively). The 
large-scale toroidal field is therefore mainly organized as an 
axi-symmetric ring, displaying a striking similarity with the 
"sea-snake" toroidal field structures produced by 3D MHD 
simulations in the presence of rapid rotation (Brown et al., 
2007b). 

In the Sun, the geometrical complexity of the large- 
scale field is varying during the solar cycle (e.g., Sanderson 
et al. 2003). An axisymmetric dipole is observed close to 
the solar minimum, but a high quadrupole/dipole ratio is 
observed at solar maximum. We can naturally expect time 
variability to occur for our stellar sample as well, implying 
that the energy splitting between low-order field components 
should be monitored over several years before discussing the 
possible informations it can disclose about the underlying 
dynamo. 

4.4 Differential rotation 

One of our target (HD 190771) possesses a measurable level 
of differential rotation. The absence of positive measure- 
ments on other targets reflects the combined effect of higher 
relative levels of noise and sparser phase sampling. With a 
level of rotational shear about 2.3 times that of the Sun, 
HD 190771 is following nicely the mass-trend derived obser- 
vationally by Barnes et al. (2005) for other rapidly-rotating 
stars (note that, probably because of this observational bias 
towards high rotation rates, the solar differential rotation 
parameters fall outside the mass dependence derived by 
Barnes et al.). The same authors also report an increase of 
dQ with Q eq (despite considerable scatter in available differ- 
ential rotation measurements). Our own analysis is in global 
agreement with previous estimates derived for similar rota- 
tion rates by Doppler imaging (Barnes et al. 2005), broad- 
band photometry (Henry et al. 1995) and chromospheric 
monitoring (Donahue et al. 1996). 

Recent 3D numerical simulations investigate global 
flows in the presence of rapid rotation (Brown et al. 2007a). 
Their model simulate the convection up to 0.95 Rq and 
ignore the very top layers of the convective zone. This ra- 
dial limit could be problematic when comparing numerical 
simulations to observations, since the rotation can possibly 
undergo steep gradients just below the surface of rotating 
stars, as it does in the Sun (Schou et al. 1998). Whatever 
the conclusions of numerical simulations, we can therefore 
expect some difference with a surface estimate of differen- 
tial rotation, as derived from photospheric magnetic maps. 
Despite this possible issue, Brown et al. (2007a) come to 
conclusions in good agreement with our observations, with 
a level of differential rotation dQ of 1 Mq stars increasing 
with the spin rate by a factor of ~3 when increasing the 
rotation from 1 to 3 Q@. 

Our isolated measurement of differential rotation 
clearly deserves to be repeated on other stars, to confirm 
the shear level we observe at 3Q© for HD 190771 and derive 
similar measurements for different rotation rates. From our 
analysis, we can already propose that an increase of dQ from 
slow to fast rotation may contribute to make the Q effect of 



the stellar dynamo more efficient and may, as a consequence, 
help to build a global-scale toroidal field throughout the con- 
vective zone. 

4.5 Long-term magnetic variability 

The various magnetic properties we derive here may be 
subject to long-term variability. The monitoring of chromo- 
spheric emission reported by Baliunas et al. (1995) and Hall 
et al. (2007a, 2007b), actually shows that two of our targets 
(HD 76151 and HD 146233) follow activity cycles. 

With a period of 7 years (Hall et al. 2007b), the cycle 
length of HD 146233 is shorter than the solar cycle. The 
chromospheric activity of this solar twin was reported to 
sharply grow between 2004 and 2006 (after a regular de- 
crease between 2000 and 2004), so that our spectropolari- 
metric observations of the summer of 2007 are likely to be 
representative of a high activity state. The predominantly 
quadrupolar magnetic geometry we reconstruct at that time 
is reminiscent of the magnetic topology of the Sun at so- 
lar maximum. However, despite the fact that HD 146233 is 
known as the best solar twin among bright stars, its short 
cycle length indicates that any comparison with the solar 
magnetic topology should be considered with caution. 

The long-term chromospheric variability of HD 76151 
(Baliunas et al., 1995) suggests the existence of a 2.5 yr 
cycle (of weak amplitude compared to chromospheric solar 
variability), superimposed on a possibly much longer cycle, 
of higher amplitude and unknown period. Its last activity 
peak was recorded in 2006 (Hall et al. 2007a), so that our 
data set was collected during the following decrease of ac- 
tivity. 

No such long-term chromospheric follow-up is available 
for HD 190771 and HD 73350. The comparison with other 
stars of similar activity levels and spectral types (Baliunas et 
al. 1995) reveal however that the cyclicity undergone by their 
Call emission is generally less pronounced than in the Sun, 
with a range of possible cycle lengths (depending on the star) 
and the co-existence of different time-scales of variability. 



5 CONCLUSIONS 

From a set of stellar spectropolarimetric observations, we 
report the detection of surface magnetic fields in a sample 
of four solar- type stars. Assuming that the observed vari- 
ability of the polarimetric signal is controlled by stellar ro- 
tation, we establish the rotation periods of our targets, with 
values ranging from 8.8 d (for HD 190771) to 22.7 d (for 
HD 146233). One star of our sample (HD 190771) displays a 
measurable level of surface differential rotation, with a rota- 
tional shear about 2.3 times stronger than solar. Apart from 
rotation, fundamental parameters of the selected objects are 
very close to the Sun's, making this sample a practical basis 
to investigate the specific impact of rotation on magnetic 
properties of Sun-like stars. 

We reconstruct the large-scale magnetic geometry of the 
targets as a low-order (£ < 10) spherical harmonics expan- 
sion of the surface magnetic field. From the set of magnetic 
maps, we draw two main conclusions, (a) The magnetic en- 
ergy of the large-scale field increases with rotation rate. The 
increase of chromospheric emission with the mean magnetic 
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Table 3. Magnetic quantities derived from the set of magnetic maps. We list the mean unsigned magnetic field (B mean ), the fraction 
of the large-scale magnetic energy reconstructed in the poloidal field component and the fraction of the poloidal magnetic energy in the 
dipolar (£ = 1), quadrupolar (I = 2) and octopolar (I = 3) components. In the last column, we also list logi?^ K values derived from our 
sets of Stokes I spectra. 



Name 




pol 


. en. 


dipole 


quad. 


oct. 


1°§ ^hk 




(G) 


(% 


tot) 


(% pol) 


(% pol) 


(% pol) 




HD 146233 


3.6 ±1 


99.3 


±0.2 


34 ±6 


56 ±6 


10 ± 10 


— 4.85 ± 0.02 


HD 76151 


5.6 ±2 


93 


±6 


79 ±13 


18 ±8 


3±3 


— 4.69 ± 0.02 


HD 73350 


42 ±7 


52 


±3 


24 ±5 


29 ±8 


33 ±5 


-4.48 ±0.02 


HD 190771 


51 ±6 


34 


± 1 


43 ±8 


20 ±2 


23 ±4 


-4.42 ± 0.02 



field is flatter than observed in the Sun. Even if we cannot 
completely rule out the possible influence of increased spot- 
tedness in more active stars, this observation may suggest 
that a larger fraction of the surface magnetic energy is stored 
in large scales as rotation increases (since the chromospheric 
flux is also sensitive to magnetic elements smaller than those 
contributing to the polarimetric signal), (b) Whereas the 
magnetic field is mostly poloidal for low rotation rates, more 
rapid rotators host a large-scale toroidal component in their 
surface field. From our observations, we infer that a rota- 
tion period lower than ^12 days is necessary for the toroidal 
magnetic energy to dominate over the poloidal component. 

In future work, we will investigate how this rotation 
threshold may vary for various stellar masses, and therefore 
study how the depth of the convection zone may influence 
the existence of global-scale toroidal magnetic fields at the 
stellar surface. This enlargement of the available stellar sam- 
ple is already under way. The magnetic properties we derive 
here may also be subject to some level of time variability. 
This idea is supported by the long-term monitoring of chro- 
mospheric emission of Baliunas et al. (1995) and Hall et al. 
(2007a), in which two of our targets (HD 76151 and HD 
146233) are reported to follow activity cycles. Variabilities 
of the large-scale magnetic geometries related to such fluctu- 
ating activity remain to be unveiled. A long-term observing 
program is now started to gather the spectropolarimetric 
data sets necessary to tackle this question, with the objec- 
tive to offer a completely new set of observables with which 
to constrain stellar dynamo theories. 
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